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^ ■ Abstract 
m 

Using the 0(4) linear a model, we address the topic of non-equilibrium relax- 
| ation of an inhomogeneous initial configuration due to quantum and thermal 

fluctuations. The space-time evolution of an inhomogeneous fluctuation of 
the condensate in the isoscalar channel decaying via the emission of pions in 
the medium is studied within the context of disoriented chiral condensates. 
Ph! We use out of equilibrium closed time path methods in field theory combined 

r~| ! with the amplitude expansion. We give explicit expressions for the asymp- 

totic space-time evolution of an initial inhomogeneous configuration including 
the contribution of thresholds at zero and non-zero temperature. At non-zero 
temperature we find new relaxational processes due to thermal cuts that have 
no counterpart in the homogeneous case. Within the one-loop approximation, 
we find that the space time evolution of such inhomogeneous configuration 
out of equilibrium is effectively described in terms of a rapidity dependent 
temperature T(i?) = Tj cosh[#] as well as a rapidity dependent decay rate 
V(d,T(&)\ This rate is to be interpreted as the production minus absorp- 
tion rate of pions in the medium and approaches the zero temperature value 
at large rapidities. An initial configuration localized on a bounded region 
spreads and decays in spherical waves with slower relaxational dynamics at 
large rapidity. 
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I. INTRODUCTION 



The dynamics of relaxation of inhomogeneous field configurations out of equilibrium is 
an important problem and a common theme in cosmology and high energy physics. 

In high energy physics the experimental possibility of studying the chiral and quark- 
gluon phase transition with high luminosity hadron colliders and upcoming heavy-ion col- 
liders makes imperative the understanding of relaxation and transport processes in extreme 
environments 

A very exciting possibility |3|-[7| that can be studied in high energy-high luminosity 
hadron collisions or relativistic heavy ion collisions within the energy range of the upcoming 
RHIC and LHC colliders, is that after the chiral phase transition, regions in which the chiral 
condensate is misaligned with respect to the vacuum state are formed. 

In these ultra-high energy heavy ion collisions (-y/i > 200Gev/nucleon) a large energy 
density (a few Gev/fm 3 ) is deposited in the collision region corresponding to temperatures 
above the critical value for chiral symmetry restoration « 200Mev. In this situation, it is 
possible that within a volume of a few fm 3 an inhomogeneous coherent field configuration is 
formed corresponding to a region in space in which the chiral symmetry is restored. 

The idea |3|-[7| is that as these regions cool down they relax towards the equilibrium 
situation emitting a large number of soft pions. Such a configuration has been dubbed a 
'disoriented chiral condensate' (DCC). 

Bjorken, Kowalski and Taylor proposed a 'baked alaska' scenario in which this con- 
figuration relaxes via the copious emision of pions strongly correlated in isospin. Such a 
scenario could also explain the Centauro events observed in cosmic rays in which the ratio 
of charged to neutral pions is different from 1/3 . 

A microscopic, field theoretical description of the dynamics of relaxation of these field 
configurations with typically large amplitudes faces at least two major obstacles. The first is 
a non-perturbative treatment of inhomogeneous, large amplitude field configurations. The 
second is a consistent description of non-equilibrium processes that can allow a real time 
calculation of the evolution. 

Recently there has been a surge of interest in the description of non-equilibrium processes 
both in cosmology and high energy physics. We refer to non-equilibrium processes to those 
corresponding to the the time evolution of quantum states in field theory which are not 
hamiltonian eigenstates as well as the time evolution density matrices that do not commute 
with the Hamiltonian. In particular the relaxation of strong electric fields H and the 
evolution of homogeneous scalar order parameters during phase transitions in Minkowski 
and cosmological spacetimes |J. 

The non-equilibrium aspects of disoriented chiral condensates have been studied within 
different approximation schemes by several authors |10|-|T^], but mainly focusing on homo- 
geneous expectation values or correlation functions. 

However, to the best of our knowledge, these studies focused on the evolution of a 
homogeneous but time (or proper time) dependent mean field and fluctuations around it but 
did not address the relaxational dynamics of inhomogeneous field configurations including 
quantum and thermal effects. 

Moreover, inhomogeneous field configurations appear in cosmology when topological ob- 
jects, such as textures or cosmic strings involving inhomogeneous field configurations are 
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considered. Their relaxational dynamics is thought to have a bearing on the spectrum of 
fluctuations in the cosmic microwave background radiation [17,18 . 



The focus of this article is to study the linear relaxational dynamics of inhomogeneous, 
coherent field configurations both at zero and finite temperature in the 0(4) linear sigma 
model. As is known such model is the relevant field theory for the description of chiral 
condensates. 

We use out of equilibrium closed time path methods flHH^B combined with the amplitude 
expansion [plj25|| to study the time evolution of (<$>(a!,t)). We consider small initial field 



amplitudes such that we can keep just the linear term in the amplitude expansion. In such 
approximation the field evolution equations linearize. 

Whereas the focus is on the description of the dynamics within the setting of the chiral 
phase transition, the analysis presented below applies to a wide variety of physical situations 
in which a inhomogeneous scalar order parameter relaxes to the equilibrium situation via 
the production of lighter scalar fields in a medium. 

The main idea is that after a phase transition (within the framework of DCC's the 
chiral phase transition) regions are formed within which the scalar order parameter is inho- 
mogeneous. As the zero momentum component rolls towards the equilibrium value, these 
inhomogeneous configurations will relax in such a way that the spatial gradients of the 
field configuration become smaller so as to decrease the energy. This relaxation will be 
accompanied by the production of quanta, such as pions in the case of DCC's. 

Here we study the 0(4) linear sigma model in the broken symmetry state, at temperatures 
below the chiral phase transition as an effective low energy theory of S77(2)l x S77(2)r (up 
and down quarks) which presumably incorporates the effects of strongly interacting QCD 
on chiral dynamics |TOf . 

Wilczek and Rajagopal PJIO| have argued that the 0(4) linear sigma model effectively 
describes the same equilibrium universality class of QCD with two flavors of light quarks 
and proposed to study it within the context of DCC. Clearly this simple model misses a 
great number of degrees of freedom, hadrons, vector mesons, etc. and can only be justified 
as an effective description at temperatures well below the chiral phase transition when the 
degrees of freedom with much higher masses are not relevant. 

The linear sigma model may also be obtained as a Landau- Ginsburg effective theory from 
a Nambu-Jona-Lasinio model PB which is a popular description of the phenomenology of 
chiral symmetry at the quark level and which has also been used as an alternative description 
of DCC's |27|. 

The 0(4) linear sigma model has also been studied within the context of texture type 
configurations deemed relevant in certain cosmological models of structure formation | |17| . 
hence our study may be relevant in these cosmological scenarios as well. 

We will restrict our study to the case in which the condensate occurs in the isoscalar 
channel, without isospin violation. In section II we introduce the model, discuss its range of 
validity and summarize aspects of non-equilibrium field theory relevant to the calculation. 
Section III presents the bulk of the results, with explicit analysis of the space-time evolution, 
including the contribution from thresholds and thermal cuts. We summarize our results in 
the conclusions section. 
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II. THE MODEL AND THE TECHNIQUES 



We study is the 0(4) linear sigma model as an effective low energy description of chiral 
symmetry aspects of QCD with two (light) flavors of quarks. The Lagrangian density is 
given by 

C = l -d$ ■ - A($ • $ - flf + ha , (2.1) 

where $ is an 0(4) vector, $ = (er, ff). The field a describes the < qq > condensate while 
the (isospin) triplet 7? describes the pions. The small magnetic field h ~ (120Mev) 3 accounts 
for the explicit symmetry breaking arising from the up and down quark masses and gives the 
pions a small mass m v pa 130Mev. The pion decay constant f n ~ 90Mev and A is fixed by the 
value of the 'mass' of the a field M a pa 600Mev pa 2v^2A/ 7r . In what follows we will neglect the 
magnetic field, but explicitly give a mass to the pions and use f n as the pion decay constant. 
At this point we recognize that the linear cx-model is a strongly coupled effective theory 
and any perturbative approximation will be subject to criticism. This important objection 
notwithstanding, we pursue the study of the non-equilibrium evolution of the condensate 
with the hope of learning new features of the decay of the condensate that will perhaps be 
rather generic and persist in higher order calculations or eventually in a non-perturbative 
treatment. 

Since we want to understand the relaxation of inhomogeneous perturbations of the con- 
densate around the equilibrium value f w , we write 

a(x,t) = f n + (f>(x,t)+x(x,t), (2.2) 
<j>(x,t) = {(a{x,t) - U)) , (2.3) 
( x (x,t)) = 0, (2.4) 

where the average above stands for the expectation value in the non-equilibrium state dis- 
cussed below, and the fluctuation field x describes the a mesons in the broken symmetry 
state. The quantum state that leads to (|2~3l) is an inhomogeneous coherent state of a mesons. 



A. Out of Equilibrium Techniques 

The field theoretical methods to describe processes out of equilibrium are known and 



described at length in the literature [|X9| [23[1 . The basic ingredient is the time evolution of an 
initially prepared density matrix, which leads to the generating functional of non-equilibrium 
Green's functions in terms of a path integral representation along a complex contour in time. 
The contour involves a forward time branch, a backward time branch and a third branch 
down the imaginary time axis to time r = —i(3 if the initial density matrix describes an 
ensemble at initial temperature 1/(3. It can be proven (see references above) that the third, 
imaginary branch only determines the boundary conditions on the Green's functions, but 
does not enter in the calculation of real-time correlation functions. 

The fields living on the forward and backward branches will be labelled with + and — , 
respectively, and the effective Lagrangian that enters in the path integral representation of 
the non-equilibrium generating functional is given by 
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£ noneq = £[$+] -C[$~]. (2.5) 

From this path integral representation it is possible to construct a perturbative expansion 
of the non-equilibrium Green's functions in terms of modified Feynman rules: 

1. The number of vertices is doubled. Those in which all the fields are on the + branch 
are the usual interaction vertices, while the vertices in which the fields are on the — 
branch have the opposite sign. 

2. The combinatoric factors are the same as in usual field theory. 

3. The spatial Fourier transform of the (bosonic) propagators are 



Gt + (t,t') = G>(t,t')Q(t-t') + G<(t,t')Q(t'-t) , (2.6) 

cr(M') = G>(t,t')e(t' -t) + Gt(t,t')e(t - 1) , (2.7) 

Gt-(t,t') = -G<(t,t'), (2.8) 

G k + (t,t') = -G>(t,t'), (2.9) 

G>(M') = ijd 3 x ($(f,t)$(0,t')> , (2.10) 

G%(t,t') = ijd 3 x e-^ ($(0, t)) , (2.11) 



where $ is a generic bose field. That is $ = cr or $ = tt. 

Now we have to specify the properties of the initial state. A particularly convenient 
choice is that of a thermal initial state at temperature T. Then the density matrix of this 
initial state is p = e~ H( ^ T , where H is the Hamiltonian for times t < 0. This choice of 
the initial state determines the boundary conditions on the Green's functions; these are the 
usual periodicity conditions in imaginary time (KMS conditions): 

G< (x, t; x 1 , t 1 ) = G> (x , t-i/3;x',t') . (2.12) 

Finally, the free-field Green's functions for generic bose fields are constructed from the 
following ingredients: 

G^(M') = ^~ {[1 + n b (u k )}e-^-^ + n 6 (u, fe )e^-*'n , 

Gt(t,t') = -M[l +n b (uj k )}e^ t ~^ +n b {u k )e~ l ^t-t>)\ 
Zuj k >■ J 

u k = \jk 2 + m 2 , n b (uj k ) = - , (2.13) 

where m is the mass of the boson. An important property that will be used in the calculations 
that follow is the relation 

G>(t,t')=G<(t',t). (2.14) 

Our goal is to study the dynamics of the expectation value of the chiral order parameter 
4>(x, t) = ((cr(x, t) — /tt)). We obtain the equation of motion for using the tadpole method 
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To implement this method in the non-equilibrium formulation, set 

a ± (x,t) = f 7T + <f){x,t) + x ± {x,t) (2.15) 



in the Lagrangian |2.5| , and consider (f>(x, t) as a c-number background field, keeping the 
linear, cubic and quartic terms as interactions. The expectation value of the fluctuation 
X ± is computed in presence of the background of fields <p and the condition {x^) = is 
imposed to all orders in perturbation theory. We restrict our study to the case of linear 
relaxation, so that we obtain the equation of motion linearized in the amplitude of (p(x,t), 
an approximation that is valid for small amplitude fluctuations from the minimum of the tree 
level potential. This equation can be obtained in a systematic loop expansion using the usual 
Feynman rules with the doubled interaction vertices and the non-equilibrium propagators. 
To one loop order there are two tadpole contributions to the equation of motion: one 
from the quartic vertex that is absorbed in a mass renormalization and the other from the 
cubic vertex that is absorbed in a renormalization of /„-. The remaining one loop Feynman 
diagrams that contribute to the equation of motion are shown in Figure 1. The details have 
already been presented elsewhere P^| , p5|| . The equation (x + (x, t)) = leads to 



d 3 x'dt'{(( X + (x,t)x + (x,t')) - ( X + (x,t) X -(x,t'))) (2.16) 
4>(x* , t') - V 2 0(f , t') + M 2 (T)0(f , t') + f d 3 x"dt"Z ret {x' - x", i! - t")<j){x", t") 



= . 

We have absorbed a momentum independent but temperature dependent tadpole in a 
renormalization of M a , and used the property ( |2.14 ). S ret is the retarded self-energy. 



We are interested in the decay of an inhomogeneous chiral condensate into pions. In the 
linear a model the chiral condensate is represented by the a field. The dynamics of the decay 
process of such configuration will be determined by the imaginary part of the self-energy on 
mass shell. 

In a medium, different kind of processes contribute to the imaginary part of the on-shell 
self-energy of the external particle |2^-^TJ. In particular, collisional processes are always 



present and are responsible for a collisional lifetime of the particle in the medium. Decay 
(and recombination) processes are only present if the kinematics of decay is allowed; these 
contribute to an imaginary part of the self-energy on shell only when the lowest multiparticle 
threshold is below the single particle particle pole in the spectral density of the one-particle 
propagator. Since the 'mass' of the a field (~ 600 Mev) is larger than twice the pion mass 
(^ 130 Mev), the one-loop diagram with pions in the internal loop will contribute to an 
imaginary part on mass shell of the sigma particle through processes in which the a decays 
into two pions in the medium, as well as the inverse process of pion recombination into o. 

Decay and recombination processes will also be present in higher order contributions to 
the self-energy, but we will only consider the contribution from the one-loop diagrams shown 
in figure 1. 

The bilinear and trilinear interaction vertices needed for the one loop calculation are 
obtained from the following part of the interaction Lagrangian density 

A = 9 {<P(x + f + (X + ) 3 + 0(v? + ) 2 + X + (??+) 2 -(+--)}, (2-17) 
M 2 

g = 4\U = -f. (2.18) 

AJtx 
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The retarded self-energy to this order is found to be: 

^>ret(x X 1 } t t ) = S. re t j( j(x x' , t t ) 4~ S re t j7r (x x' , t t ) , 

E^^f-x'^-t') = 18 ig 2 {[G>(x- x',t-t')] 2 - [G<(x - x' ,t - t')f}Q(t - t') , 
E ret ^(x-x',t-t') = 6 ig 2 {[G>(x-x',t-t')} 2 - [G<(x - x' , t - t')] 2 }Q(t - t') . (2.19) 

The Green's functions in ( |2.19| , |2.19|) are given by ( |2.13|) appropriate for the quanta of 
the fields \ and 7? respectively. We have also used the fact that the pion Green's functions 
are diagonal in isospin. 

Since we are dealing with a real field, the retarded self-energy E ret is given by [ [22] , f2"3[l 

E ret (£ , t;x',t') = 2 ReY? (x, t; x*, t') . (2.20) 

Translational invariance of the self-energy makes it convenient to write the equation of 
motion for the spatial Fourier transform of 4>. Introduce 

mt) = /t|^ c ** ( 2 - 21 ) 

E ret (f - x , t - t') = I eW*-*) E(p, t-t'). (2.22) 

In order to solve the integro-differential equation we will impose the initial condition 
(fi(x, t < 0) = and that this configuration is 'released' at time t = f32fl . Under these 
conditions the evolution equation becomes 

8(p, t) + u 2 p 5{p, t) + J dt' E(p, t - t') 5{p, t') = 0, (2.23) 
u 2 p = \p\ 2 + M 2 (T) . (2.24) 

Using the free-field Green's functions given above for the x > ^ fields we find 

E(p, t-t') = E a (p, t-t') + E ff (p, t - t') 

r d 3 k 1 

T, a (p,t - t') = -18g 2 / — r-r — {(l + 2nk,a)sm[(uk+p, <7 + ^k,a)(t-t , )\ 

-2 n Ku sm[(u k+P:a - u Ka ){t - t'))} (2.25) 
d 3 k 1 
(27r) 3 2u kj7r aj k+Pi7r 

2 n k)lv sm[(u k+Pt7T - u k ^) (t - t')) } (2.26) 



r d 3 k 1 

E w (p,t - t') = -Qg 2 / — — -- {(1 + 2n ky7T ) sin[(uj k+p ^ + w k v )(t - t')] 



where 



u k ,w = sjk 2 + M CT % , (2.27) 

*W = e **.l - 1 ' ^ 2 - 28 ) 

The equation of motion (|2.23|) can now be solved via a Laplace transform. With the 
boundary conditions S(p, t — 0) = 5i(p) ; <5(p, £ = 0) = and denoting the Laplace transforms 
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of S(p,t) , E(p, t) by 5(p; s) , S(p; s) respectively (here s is the Laplace transform variable) 
we find 

^»= ?4fW (2 - 29) 

with u;^ given by ( |2.24| ). The time evolution is found by integration in the complex s-plane 
along the Bromwich contour s = iu + e ; —oo < u < oo with e chosen so that the contour 
lies to the right of the real part of all the singularities of 5(p; s). 
The self energies (|2.25| , |2~2~6"1) are proportional to the kernel 

( d 3 k 1 

K(p, t-t') = -2 {(1 + 2n k ) sin[(o;^+ u> s ){t - t')] 

-2n k sm[(u; n+f -^)(t-t')}} , (2.30) 
whose Laplace transform can be written as a dispersion integral 



CO 



K(p;s) = -I dVo ^TV Y ' ( 2 - 31 ' 



2 p p(p,p ) 
s 2 +pl 
[ d 3 k 1 

p(p,Po) = / / 9 x 39 , — {(! + 2n fc )5( Po - ^ - - 2rz fc 5(p - ^ + ^ +p .)} . (2.32) 



(2vr) 3 

In the expressions above, the frequencies and Bose factors correspond to the particle in the 
loop. 

We are considering the situation near the end of the chiral phase transition, in the 
hadronization stage with temperature different from zero but well below the transition tem- 
perature < T < M v . However, the zero temperature limit can be taken at any stage in 
the calculation. 

Since this kernel will determine the singularities of the Laplace transform in the com- 
plex s-plane it is important to understand its analytic structure. From (|2.32|) we find the 
imaginary part of the kernel to be 

ImK{p,iuj ± + ) = ±7rsign(co>) [p(p, — p(p, — , (2.33) 

where the sign(u;) reflects the retarded nature of the kernel P3LI30 . 



There are two different processes that contribute to the imaginary part of the kernel that 
deserve to be studied in detail, for u > these give the following contribution: (for uj < 
changes sign) 

ImK(p, ioo + + ) = ImK {1) (p, iu + + ) + ImK {2) (p, iu + + ) , (2.34) 

1 ( d 3 k 

ImK^ip, iui + + ) = — 2 (1 + 2n k )S{u - u n - , (2.35) 



ImK {2) (p, iu + 0+) = -— / 2n k 5(u + u u - u % .) -6(u e -u- . (2.36) 

107T J UJ^UJ^ +p 



The processes that contribute to ImK^ p9|-pT| are the decay a — > nix with Boltzmann 
weight (1 + n£)(l + n^ + ^) minus the recombination process tttt — > o with Boltzmann weight 



S 



n k n k+p wnen the particles in the loop are pions, or the process a — > 2a and its inverse 
2a — > a with similar weights when the particles in the loop are a. Clearly these processes 
will only contribute above the two-particle threshold. 

The processes that contribute to ImK^ are the <77r — > 7r with weight nr(l+nt , -) 

minus the inverse it — ► an with weight (1 + n^)n^ + ^ when the particles in the loop are pions, 
or the process aa — > a minus its inverse a — ► ao" with the corresponding weights when the 
particles in the loop are the a mesons. Clearly these latter process can only occur in the heat 
bath for non-zero momentum transfer and will lead to new thresholds and discontinuities. 
These are identified with the Landau damping processes in the medium JJIJ as they only 
occur for space-like momenta. An important point to consider is that these processes do not 
contribute to the relaxation of a homogeneous condensate (zero momentum transfer) and 
are, therefore, a new feature of the inhomogeneous situation. 

After analyzing the kinematical regions, we find: 



ImK {1 \p,iuj + + ,T) -- 
x6(w 2 -p 2 - 4m 2 
ImK {2 \p,iuj + + ,T) = 



1 



57T 



1- 



Am 2 



UJ 



P A 



T , 

+ - — In 

Anp 



I — e -/3u p 



T 

Anp 



In 



®{p 2 




4m 2 



uj 2 — p 2 



(2.37) 
(2.38) 

(2.39) 



where m is the mass of the particles in the loop, m = M v for or m = M a for S CT . In 
the expression for ImK^ 1 ' we have explicitly separated out the T = and T / contri- 
butions. Note that whereas IraK^ 1 ' is non-zero above the two particle threshold and has a 
contribution that survives in the zero temperature limit. ImK^ is non-zero only for 
and for space-like momenta and does not contribute directly to the imaginary part on shell 
but does contribute to the relaxational dynamics as discussed below. The explicit form of 
the imaginary part at finite temperature is one of the novel results of this study. For small 
momentum p/T < 1 we find 



ImK {1) {p,iu + + ,T) 
ImK ( - 2 \p,ito + + ,T) 



ImK (1) (p, iui + + , T = 0) [l + 2n(u; 




Am 2 



;2n{uo-)Q{p 2 - uj 2 ) 



(2.40) 
(2.41) 



and we recognize (|2.40| ) as the imaginary part obtained in the homogeneous case |25 |. 
The 7r ; a self-energies are now given in terms of the above kernel as: 

S ff (p, uj, T) = 3g 2 K{p, uj;m = M n ,T) , (2.42) 
E CT (p, uj, T) = 9g 2 K(p, uj;m = M a ,T). (2.43) 
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B. Zero Temperature Limit 



The real part of the kernel can be obtained from the imaginary part from the dispersion 
integral ( |2.31| ), but is very difficult to compute in the general case. Its zero temperature 
limit can be found analytically. The integral in ( [2.31|) has a logarithmic divergence which is 
subtracted at s = ;p = 0. This subtraction can be absorbed in a further renormalization 
of M a and we find that for T = 0: 



K(p,s) 



An 2 



'1 



Am 2 
s 2 + p L 



Argtanh 



1 



/ -r _ 4m 2 
1 ~r «2 , „2 



1 



[2.AA) 



where again, m refers to the mass of the particle in the loop. Along the imaginary axis and 
below the two particle cut, the kernel is real and given by 



K(p,s 



ILO 



4vr 2 



Am 2 



— 1 arctan 




- 1 



(2.45) 



whereas above the two-particle cut uo 2 > Am 2 + \p\ 2 the kernel has an imaginary part, 
K(p, iuj ± 0+) = K R (p, u) ± iKiip, u) with 



Kr(p,uj) 



An 2 



Am 2 
uj 2 — p 2 



Argtanh 



4m 2 
uj 2 — p 2 



- 1 



1 



4m 2 
uj 2 — p 2 



Using the parameters of the linear a model, we find that 
ReE^(p,u,T= 0) +ReX a {p,uj,T = 0)| 
ImTiK^p, u, T = 



1.195Mf . 



0.045il€ 



(2.46) 

(2.47) 
(2.48) 

(2.49) 
(2.50) 



Thus we see that although the linear a model is strongly coupled, the one-loop correction to 
M a is very small. However the imaginary part on shell (related to the width of the particle, 
see below) is very large. 



C. Analytic structure 

To perform the inverse Laplace transform we need to understand the analytic structure 
of the Laplace transform in the complex s-plane. 

From the analysis presented above, we find that for arbitrary temperature, 5(p; s) is 
analytic except for 

1. isolated single particle poles at s — ±iQ(p,T) with Q(p,T) being the solutions of 

- Q 2 (p,T) + \p\ 2 + Ml{T) + E(s = ±ifl,p,T) = (2.51) 

and 
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2. discontinuities across the imaginary axis determined by the imaginary part of the 
self-energy obtained above: 

5(p, s = iu + + )- 5(p, s = ito-0 + ) = i 5i(p) S(u,p,T) , (2.52) 

where we have introduced the spectral density S(u,p,T) given by 

^ 2u £/(u;,p, T) 

S[UJ ^ T)= p - \pf - Ml - Y^pVrf + ^,P,T) 2 ' (2 ' 53) 



III. TIME EVOLUTION OF THE ORDER PARAMETER 

The real-time evolution of the spatial Fourier transform of the condensate is obtained 
via the inverse Laplace transform 

S (0 t ) = — e st S Ml (3 1) 

The integral is performed by deformation of the contour, wrapping around the single 
particle poles which are solutions of the equation ( |2.51|) and slightly to the right and left 



of the multiparticle cuts along the imaginary axis. Thus the inverse Laplace transform will 
have two contributions: from the poles S po i e (p,t) and from the cuts S cut (p,t) and we can 
write in general 

5{p, t) = 5 po ie(p, t) + 5 mt (p, t) ■ (3.2) 
There are several different cases that we can study. 



A. Stable Case (M a < 2M W ) 

Although this is not a realistic case for the chiral phase transition described by the 
linear a model, this case is still worth studying in order to compare it to the unstable 
case [M a > 2M 7r ) to be described later. Furthermore, we will see that many results from 
this particular case do apply in the unstable case which is of interest for the chiral phase 
transition. 



1. Zero Temperature 



At zero temperature the self-energy is manifestly Lorentz covariant as can be seen from 
the expressions ( [2.46 , |2.47p . 

We first analyze the one-particle pole contribution, which is a straightforward general- 
ization of the homogeneous p = case [p5 |. The condition (|2.51|) is satisfied for 



fi 2 (p,T = 0) = |pf + M 2 , 
M 2 -M 2 -S(s = iM ,p = 0,T = 0) = 0. 



(3.3) 
(3.4) 
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Then the pole contribution to the inverse Laplace transform is 



S P oie(p, t) = 5i(p)Z cos( yV + Ml t) 
dZ(iM,p = 0) 



dM 2 



M=M„ 



(3.5) 
(3.6) 



where Z is the (finite) wave-function renormalization constant, defined on shell. We will 
now consider the case in which Si(p) = 5, (independent of p) and later study the general 
case by convolution. 

At this point it is convenient to introduce the proper time and (radial) spatial rapidity 
variables as these exhibit the Lorentz properties more clearly: 



iln 

2 



t + r 



t — r 



T 



rsinhf^] , t = Tcosh[-# r 



The spatial Fourier transform leads to the pole contribution 

Z5i Ml cosh[$ r 



</W(^, t) 



2tt 



T 



J 2 (M t)9(t 2 



(3.7) 



(3.8) 



where J 2 is a Bessel function. Recalling that bi = J d 3 x<j)(x, t — 0) it is clear that this solution 
is Lorentz invariant, since the product 5i cosh[$ r ] is Lorentz invariant. The pole contribution 
then gives fluctuations which propagate inside the light-cone as a massive relativistic wave 
as expected. For large times and distances eq.(|3.8[) gives 



4>pole(x,t) 



r.t— >oo 



Z5, 



-M, 



3/2 



cosh[-^ r 



cos(M, 



(3.9) 



^3/2 « r 3/2 — v-w 4 . 

We now evaluate the cut contribution 6 cut (p,t) by performing the integral over the cut and 
using the change of variables uj —>■ \/uj 2 + p 2 



8i r 00 1 

o~cut(p,t) = — du S(u,p= Q,T = 0) cos(Juj 2 +p 2 t) 

7T J2M 



(3.10) 



where M. is the smallest of M a ; M n and S(u,p,T) is the spectral density given by ( p.53|) . 
The amplitude is then 



4>cut(x,t) 



5i cosh [i? r 



4tt 2 



/ u 2 S{u,p = 0,T = 0) J 2 {ujt) du . (3.11) 

J2M 



For large proper time r ^> M n 1 ^> M a 1 this integral is dominated by the two-particle 
thresholds (uj ~ 2M CT ; uj ~ 2M n ) and we find 



, r,t->oo , COSh[l? r ] 
4>cut{x,t) = bi 5 



CiM^ sin(2M (J r) + C 2 M 2 sin(2M 7r r) 



(3.12) 



where C\ ; C2 are constants determined by the spectral densities at the respective two-particle 
thresholds. When Ml r 2 ; M 2 r 2 ^> 1, the amplitude of 4> cut is smaller than the one of 4> po i e 
by a factor r 3 / 2 and is clearly subleading at long distances. This is analogous to the factor 
t 3 / 2 for the homogeneous case p5 . 
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2. Non-Zero Temperature 



Unlike the T = case, we cannot make use of manifest Lorentz covariance, since this 
property is not available in the rest frame of the heat bath. However the situation is 
qualitatively similar to the zero temperature case. The positions of the one-particle pole are 
now given by 



s = ±iQ(p, T) 



Sl(p, T) = y Ml + p 2 + £(ifi(p, T),p, T) . 
The pole contribution is 

SpoUf, t, T) = SiZ(p, T) cos (fi(p, T) t) , 
where the wave-function renormalization (defined on shell) is now 

dE(iM,p,T) 



Z(P,T) 



dM 2 



M=n( P ,T) 



(3.13) 
(3.14) 

(3.15) 
(3.16) 



The lack of Lorentz invariance introduces a complicated dependence on the momentum p, 
and the inverse Fourier transform cannot be performed analytically as it was the case for 
T = 0. 

However, for t > r ^> M" 1 , M~ l (t > r 3> O(lfm)) one can compute (ppoiei^, t, T) by the 
stationary phase approximation. We have 



4>pole{x,t,T) 



Si 



2n 2 r Jo 

This integral has stationary points at 



dpp sin(pr) Z(p, T) cos [t T)} . 



M„r 



y/t 2 -r 2 



+ 0(g* 



(3.17) 



(3.18) 



for large t and r and time-like intervals, where by 0(g 2 ) we refer to terms that are small in 
the formal weak coupling expansion or in the strongly coupled a model numerically small 



as given by equation |2.49 



These saddle point values have an illuminating physical interpretation which is better 
displayed in terms of (radial) rapidity ( |3.7|) and momentum rapidity variables defined by: 



B = - In 
p 2 



uj+p 
u — p 



p = M a sinhf^p] , uj = M a cosh^p] 
Using these variables, the saddle point condition becomes 



(3.19) 
(3.20) 

(3.21) 



and we will just refer to •& as the rapidity variable irrespective of coordinate or momentum. 
The contribution from the saddle points yields 
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Vpole 



[x, t, T) 



r,t— >oo 



6i 



X COS < M a T 



V^TT 3 / 2 
1 



Z(p = M a sinh[$],T) 



M 3 / 2 cosh[tf] 



1 + 



r 3/2 



+ 0(9* 



2M 2 
1 1 1 



E(iM,p = M CT sinh[tf],T) +0(g 



71 

4 



{ ^t 2 ' r 2 ' tr^j 



(3.22) 



The fluctuations for T ^ propagate as spherical waves, similarly to the zero temperature 
case displayed in eq. (E>.9|). However the phase and amplitude depend on r and $ through 
the temperature dependence of £ as well as of Z, which breaks manifest Lorentz covariance. 
The cut contribution is 



5 cut (p, t, T) = 8 cut ,i(p, t, T) + 8 cut>2 (p 7 t, T) , 
where we have separated the 'normal' (zero temperature) two-particle cut 

5 C ut,i{p,t,T) = — . S(u,p,T) cos(u t)du 

71 Jy/4M 2 +p 2 

from the thermal contribution with support below the light-cone 



(3.23) 



(3.24) 



S cut2 (p,t,T) 



5 { pp 



TT JO 



S(cu, p, T) cos(a; t)du 



(3.25) 



Again, as in the zero temperature case M. is the smallest mass. 

Let us first study 5cut, i- At large t the most important contribution to the integral will 
be from the region near the thresholds and the largest contribution will be from the smallest 
threshold (two pion) . From the expression for the imaginary part given by eq. fl2.37| . pT3D| ) 
we see that near the thresholds uj + — > and 



l + 2n 



(3.26) 



with n the Bose occupation number. Performing the change of variables lu — > ^cu 2 + p 2 the 

integral over to is similar to (|3.11|) but with the finite temperature factor 1 + 2n f v^L^2 \ 

in the integrand. For large distances, the spatial Fourier transform can be done by the 
stationary phase approximation. The values of the stationary phase are 

p s = ±cjsinh[tf] . (3.27) 

Finally, at large times the integral over u obtains the largest contribution near the thresholds 
to = 2M n ; uj = 2M a and we find the leading behavior for large times and distances (for time- 
like intervals) 

cuM (f , t, T) r ' =°° 6i ^M | Cl (T) [1 + 2n (M a cosh[0])] M 2 a sm(2M a r) + 

CliT) [1 + 2n (M n cosh[0])] M 2 sin(2M^r)} , (3.28) 
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with Ci(T) ;C2(T) the finite temperature counterpart of the zero temperature constants and 
the Bose enhancement factors are a result of the enhancement factors near thresholds as 
given by eq.( |3.26| ). The temperature dependence of the constants C\ ; C 2 is rather weak 
through the temperature dependence of the pole which can be neglected to the lowest order. 
Then to this order we see that the most important effects of temperature are through the 
Bose-enhancement factors which are the same as in the homogeneous case but with an 
effective, rapidity dependent temperature = Tf cosh[$]. 

The second cut gives the following contribution to the amplitude 



4>cut,2{x,t,T) 



St 



47T 3 



r Jo 



dppsmpr / S(uj,p,T) cos(ut)duj 

J—v 



(3.29) 



For large time and distances it is dominated by the point ou 
obtain 

ImZ n (pp,p,T) w "=° 3g 2 F n (p,T) + 0(p 2 ] 
ReE v (pp,p,T) "' P =° 3g 2 G T (p,T) + 0(p 2 ) 



p = 0. Defining p = -, we 



(3.30) 
(3.31) 



where the function G^p, T) can be found from the appropriate self-energy after some algebra 
but will not be relevant for our purposes, and the function F n (p, T) is given by 



Fn(p,T) 



47T 



(3.32) 



• T\/T~- 



1 



So- has an expression similar to E^, but with the factor 3 replaced by 9 and M n by M a . 

We will approximate M 2 + 3g 2 G n (p,T) by the 'pole mass' which at low temperature 
(compared to M a ) can be approximated by M a as discussed above. It is convenient to 
change integration variables in eq. (|3.29|) setting y = pr. For large r and t and T <C M a 
the contribution from the loop with a mesons is negligible, since only the lowest mass 
intermediate states contribute a long distances and the a mesons are not thermally excited. 
In this case we find 



4>cut,2{x,t,T) 



Si 3g 2 [°° 3 A ■ f +1 A 
2^M^h V dVSmy L Pdp \ 



F n (p,T) 



+ 



^F n (p,T) 



cos 



■typ. 



(3.33) 



Using the symmetry of the integrand we can recast this result as 

W r°„ 3 ^ f+\^ F n (p,T) 



2vr 3 Mi 



r jo 



p dp 



pdp 

1 1 + 



%F^T) 



sin[p(r - tp)\ . (3.34) 



The integral over p yields the third derivative of the Dirac delta and we finally obtain 



4>cut,2(x,t,T) 



r.t— *oo 



S l 3g 2 



2tt 2 M* sinh[^] cosh 4 



9{t 2 ) d 3 
dp* 



PFAP,T) 



1 + 



MS 



fAp,t) 



p=tanh [{f\ 



(3.35) 
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In the general case, there is a similar contribution from the a meson loop, but with F w 
replaced by F a as mentioned above. 

Finally the large distance, large time behavior is given by 

(f){x,t) = <f) p0 le{x,t) +<j>eut,l{x,t) + (j) C ut,2(x,t) , (3.36) 

with the asymptotic results given by ( |3.22| , E3.28"| , |3.35D . 

We emphasize that the contribution from the thermal cut yields a non-oscillatory con- 
tribution decaying as t~ 5 or r~ 5 and is subleading with respect to the contributions from 
the pole and from the 'normal' two-particle cut given by eq. (|3.28|) at long distances. Thus 



we obtain one of the important results of this work: the contribution from processes such 
as cm — ► 71 (and the inverse process), that can occur only in a plasma of excitations, are 
asymptotically subleading compared to the one-particle pole and the contributions from 
processes such as a — > ixix (and its inverse). 



3. Time evolution for an initial Gaussian wave packet 



In this subsection we will compute the time evolution in the case in which the initial 
state is a Gaussian wave packet. That is, we will assume 



[x, t = 0) 



3/2 _| 2 | 2 
I e 2 « 2 



(3.37) 



where we included a convenient normalization factor. We have to convolute this initial 
condition with the results obtained in the previous subsection. Namely 



<^(£, t,T)= / d 3 y </>t(x-y,t = 0) (f)(y, t, T) 



(3.38) 



where (p(y, t, T) is the amplitude for the evolution with the initial condition corresponding 
to a <f>(x, t = 0) = 5i5 3 (x) as studied in the previous sections. 

Within the setting of the chiral phase transition and the possibility of the formation and 
evolution of DCC's, we are interested in initial packets of spatial sizes £ > lfm (therefore 
£ 3> M~ x w 0.2fm), and the evolution for distances and times r ; t ^> £. 

We first analyze the pole contribution. By using (|3.17|) in ( |3.38|) and performing the 
simple integral over d 3 y, we get the result 



1 f°° 1 2t2 

4>t )P oie(x,t,T) = — — / dppsm(pr)e ^ * Z(p,T) cos 

Z7TT JO 



ty/p* + M* + E(in(p,t),p,T) . 

(3.39) 



We see that the only difference with (|3.17|) is the Gaussian factor in the integrand. 
For r , t ^> £ ^> M~ x the stationary phase points are 



±p = iM^sinhftf] 



1 + 



'M a e 



T 



(3.40) 
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This analysis holds for all other contributions (from both cuts), with the final conclu- 
sion that asymptotically for large distances and times compared with the size of the initial 
packet, the results obtained for the ^-function initial condition can be used with the simple 
modification of the additional Gaussian pre-factor 

Thus at time and distances much larger than £ we obtain 

<M*,t,T) = e- l ^ 2sinh2 ^ =0 (x,t,T) . (3.42) 

Only for distances comparable to the size of the initial packet or very close to the light cone 
do the corrections to the above results arising from the initial size of the packet become 
important. It is clear that the asymptotic space-time evolution is completely determined 
by the lowest energy thresholds and the fact that higher mass degrees of freedom are not 
described by the linear a model should not modify the asymptotic behavior for distances 
larger than a few fm. 

A remarkable consequence of the saddle point conditions which determine the long time- 
long distance behavior of the contributions from the multiparticle cuts is that the tem- 
perature effects enter in terms of an effective temperature that depends on rapidity given 
by 

T ef 0) = ^ rn1 (3.43) 
m ; cosh[tf] v ; 

in the sense that the Boltzmann factors that enter into the expressions have this effective 
temperature dependence. Thus the space-time evolution can be interpreted as that the inho- 
mogeneous, non-equilibrium configuration is relaxing with a rapidity dependent temperature 
that becomes smaller at large rapidities. This is a one-loop result and may be modified by 
higher order corrections, a possibility worthy of further study. 



B. Unstable Case {M a > 2M W ) 

As is well known (see [p5| , j33|1 ), if the a particle is unstable (resonance), namely if M a > 
2M n , the pole is above the two-particle threshold and moves off into the second (unphysical) 
Riemann sheet at a distance 0(g 2 ) from the imaginary axis in the s-plane. Thus in this case 
we have no pole contribution. However, if the theory is weakly coupled, the pole is very 
close to the cut, and will give the dominant contribution to the integral over the cut itself. 
In this case we have to consider only the cut contribution since, as explained in detail for 
the homogeneous case [p5|l , the poles of (pp(s) get a real part and move off into the second 



Riemann sheet in the s-plane. Thus 

<Wp,*) = 0. (3.44) 

Thus in this case we find 

6{p,t) = - 6i(p) y- — -— — — — — ^ , 12 . — ^7^- 3.45 

7T Jo [uo 2 — \p\ 2 — — Er(u,p, T) 2 + Ei(u,p, T) 2 
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This expression leads to the following sum rule obtained previously by Pisarski [34J in a 
different manner: 



" S ' ( ^' T) -1. (3.46) 



vr Jo [oo 2 - \p\ 2 - Ml - E*(u;,p,T)] 2 + E T (u,p,Ty 



In this unstable case, this sum rule illuminates the relationship between the Landau 
damping contribution to the absorptive part of the self-energy and the on-shell contribution. 
Notice that the Landau damping contribution to the absorptive part of the self-energy 
ImKW given by eq. (|2.38| ) has the same form as the finite temperature contribution to 



ImKW ( P37| ). Thus the sum rule above determines that the contribution to the spectral 
density below the light cone, that is the Landau damping term given by eq. (|2.38 ) is borrowed 



from the spectral density above the physical two-particle cut. Therefore, although ImK^ 2 ' 
does not contribute directly to the width of the meson and its decay rate in the medium, it 
does so indirectly through the sum rule. 

In weakly coupled theories such that g 2 <C 1 the spectral density features a narrow 
resonance and has the Breit-Wigner form. Thus the integral over the discontinuity feels the 
pole and can be approximated by a Breit-Wigner resonance in the second Riemann sheet. 
When the position of the resonance is far away from the two particle thresholds (many 
widths) the integral is dominated by the resonance for long times but eventually the large 
time behavior will be determined by the behavior of the spectral density at threshold [^5| . 

In what follows we will focus on the weakly coupled case and address the strongly coupled 
linear a model afterwards. 

1. Zero Temperature 



In the case g 2 <C 1, and after the change of variables uj — > \juj 2 + |p| 2 the spectral density 
can be approximated by the Breit-Wigner form 

b[LJ,p, T = 0) « , (3.47) 

{uj 2 — Mq) + Sf(a? = M ,p = 0, T = 0) 

M 2 = M 2 + E R {u = M ,p = 0,T = 0) . (3.48) 



The integral over the discontinuity of the self-energy, given by p.!0| , can be approximated 

by 



5 cut (p, t)^S,Z e~ r(p) * cos(t VX 2 + p 2 + a) , T(p) < M , (3.49) 

where the wave function renormalization Z and a p-independent and given by 



dX R (iM ,p = 0) 
DM 2 

and the damping rate is given by 



-i 



a = -Z m , (3.50) 



r&O - z ** m >p = °) . (3.5i) 



2JM 2 +p 2 
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At one loop order we find from ( [2.471 ) an d (|2.42|) 

3^ 2 



T(p) 



16nJM 2 + p 2 \ 



1 - 



4Mg _ M g r(0) 
' . / M 2 + p 2 



(3.52) 



The expression (|3.52|) displays the Lorentz contraction of the width (dilation of the 
lifetime) with respect to that in the rest frame. 

The Fourier transform of S cut (p, t) can be evaluated for large t and r with the stationary- 
phase method. We find 

5, t Z Ml cosh[tf] 



<j>(x, t) 



r.t— >oo 



-r(0)r 



x sin < M a r 



An 2 r 



+ a 



(3.53) 



There is now an exponential damping in the amplitude. However for proper times longer 
than m T _1 (0) ln[r(0)/M CT ], the most important contribution to the integral arises from the 
region near the two-pion threshold leading to a long time tail given by the pion contribution 
to 4> cut (x,t) in equation ( |3.12| ). 

For an initial Gaussian wave packet ( ^.37| ) the field <j)^(x,t) is given for large t and r by 

4>t(2,t,T) = e-^ 2si ^^ =0 (x,t,T) (3.54) 
just as in the stable case because the stationary phase condition is the same as in that case. 



2. Non-Zero Temperature 



As shown in equations ( p.23| , |3~24"| , |3.25|) , at finite temperature the are several different 
cuts which contribute to 5 cut (p,t,T). The most important contribution arises from the two- 
particle cut starting at u 2 = 4M 2 +p 2 corresponding to the process of a decay into two pions 
and pion recombination. This is clearly seen because the resonance will be at a position 
4M 2 + p 2 < Q 2 (p) < 4M^ +p 2 because the a meson can now kinematically decay 'on-shell' 
in two pions. For weak coupling the sharp resonance dominates the integral and one obtains 
a result for 5 cut) i(p,t,T) analogous to ( |3.49| ), with the finite temperature self-energy. The 
two-particle cut starting at 4M^ + p 2 and the thermal cut below the light-cone will give a 
contribution similar to the stable case studied in the previous section. Thus we concentrate 
on studying the contribution from the two-pion cut. Since manifest Lorentz covariance is 
lost (in the rest frame of the thermal bath), in the T/0 case all the parameters will depend 
on p 



Z(p,T) 



dZ R {iM,p,T) 



T(p,T) = Z(p,T) 



BM 2 
S 7 (^(p,T),p,T) 



a(p,T) = -Z(p,T) 



20(p,T) 
0S f (iM,p,r) 
dM 2 



(3.55) 

(3.56) 
(3.57) 



M=Q{p,T) 
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where fl(p,T) is given by equation ( |3.14|) . 

The inverse Fourier transform of S cutt i(p,t,T) can be evaluated for large times and dis- 
tances with the saddle point method to obtain, in the Breit-Wigner approximation 



J BW 



[x, t, T) 



5t M*Z(h(#),T) coshftf] 



-r(0,T)T 



sin[M a r + a(h,T)) 9(r) 



An 2 r 

where h(r,t) = M a sinhf^]. More explicitly for weak coupling, r($, T) takes the form 



(3.58) 



r(tf,T) 



Sg 2 T 



log 



sinh «4+(#) 



8ttM2 sinhftf] cosh[$] sinh.4_($) 



(3.59) 



with 



A±(r,t) 



M CT cosh[tf] 
4T 



1 ± tanh[tf] 



\ 



AMI 
Ml 



(3.60) 



For high temperatures and M 2 -C M 2 one finds the remarkable result (valid in the weak 
coupling case) 



r(#,T) 



3g 2 T d 



2vrM2 sinh [2^] 



(3.61) 



Notice again that the Boltzmann factors depend on the effective temperature T(i9) = 
T I cosh[#] and that T is a function oft? andT{&). In particular, notice that ( |3.59|) approaches 
the zero temperature limit at very large rapidities. 

For long proper times, r w T~ l ln(r/M CT ), the threshold will dominate the long time- 
large distance behavior with a contribution which is given by eq. Q3.28| ). The contribution 
from the thermal cut is the same as in the stable case and given by eq. (|3.35|) . In the weak 
coupling limit the contribution from the Breit-Wigner resonance and that of the threshold 
can just be added to yield the total contribution from the cut, and since the contribution of 
the thermal cut is added to that of the 'normal' cut, we finally find in the unstable case: 



4>unst(x, t, T) = (f) B w(x, t) + (j)thresh(x, t) + CJl4 , 2 (^, t) 



(3.62) 



with the asymptotic results given by Q3.58| , [3.28"l , |3.35|) 

It is now straightforward to show as in the previous section for the stable case, that for 
an initial Gaussian wave packet the final result is just multiplied by the factor e~2 M ^ smh 1*1 
because of the saddle point condition. 

Thus at time and distances much larger than £ we obtain 



'unst ;£ 



[x, t, T) 



'unst ;£=0 



% t, T) 



(3.63) 



The remarkable result of this section is that the space-time description of the evolution 
of a wave packet that decays, is that the effective decay rate is a function of rapidity which 
itself depends on the rapidity dependent effective temperature T(i9). We recall that in a 
medium this decay rate does not give the rate of production of pions through the decay of the 
sigma meson, but the difference between the production and annihilation rate as discussed 
previously in section II. 
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Let us now discuss the strongly coupled case. The approximations and conclusions 
obtained above are within the framework of a weakly coupled theory. As was discussed in 
section II, the imaginary part on shell for the linear sigma model description of pion physics 
is very large (see equation |2.50|) , thus the validity of the Breit-Wigner approximation to 
describe the space-time evolution must be called into question. Figure 2 shows S(y,p)/M a 
as a function of y — y/u 2 — p 2 /M a for T = 0.1M CT ; \p\ = 10M CT compared to the Breit- 
Wigner approximation. We see that in most of the range, but in a small region near the 
lowest threshold the Breit-Wigner approximation is excellent. We find the same result for 
temperatures all the way up to the pion mass and for a large range of momenta from very 
small to very large compared to M a . Clearly what is less accurate within the Breit-Wigner 
approximation is the extension of the integration region in u to — oo because the resonance is 
rather broad. However one can extend the integration region and subtract the contribution 
from threshold, which is the dominant one at large rapidities and proper times. This is 
precisely what was done to arrive to the result given by eq. ( |3.62| ) above. Thus we conclude 
that even in the strong coupling case the asymptotic behavior is well approximated by 
eq.(|3.63|) for the case of a Gaussian initial packet. 



IV. CONCLUSIONS 

We have studied the relaxational dynamics of an inhomogeneous condensate fluctuation 
in the 0(4) linear sigma model near the broken symmetry state, both at zero and non-zero 
temperature. 

Explicit expressions are obtained for the self-energies at zero and finite temperature and 
we point out that at finite temperature there are new relaxational processes with origin in 
thermal cuts and that are only present in the inhomogeneous case with no counterpart in 
the relaxation of an homogeneous condensate. 

For initial Gaussian fluctuations we have given explicit expressions for the asymptotic 
space-time evolution of the inhomogeneous fluctuation including the effect of thresholds both 
at zero and finite temperature. At finite temperature we obtained the decay rate of this 
inhomogeneous non-equilibrium configuration, in the medium this quantity is the rate of 
production minus the rate of absorption of pions. The space-time evolution is described in 
terms of an effective temperature and "decay rate" that depend on rapidity. We find to one 
loop that relaxational processes are described in terms of T($) = T j cosh[i9] and r($, T ($)), 
and for large rapidities the finite temperature decay rate approaches the zero temperature 
limit. 

We systematically compute the field behaviour for large times and distances compared 
with the inverse of the typical mass scale (M) in the model. This means t 3> 1/M, r 3> 1/M 
and r = \Jt 2 — r 2 ^> 1/M. For t = r ^> 1/M the behaviour will be quantitatively different. 

For large r, the field <j)(x, t) propagates as spherical waves for an initial wave packet of 
arbitrary shape concentrated around the origin. 

Our asymptotic results can be summarized as follows. 

• Stable Case: M a < 2M n 

T = 0. The pole contribution (j) po i e (x,t) dominates, giving an amplitude that decays 
as r -3 / 2 and oscillates as a function of r with frequency Mo (physical a mass). 
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The cut gives contributions smaller than 4> po i e by a factor r 3 / 2 and oscillating 
with frequencies equal to the threshold positions [eqs.( |3.9j - |3.12| )]. 

T / 0. The pole contribution dominates for large r. $ vo \e decays with the same power 
t~ 3 / 2 as for T = 0, but the amplitude and frequency change quantitatively [see 
eq.flpH)] . 

A new cut running from — ip to +ip appears at non-zero temperature. It gives a 
non-oscillating contribution for large r that decays as r~ 5 [see eq.( |3.35| ) ]. 

• Unstable Case: M a > 2M W 

T = 0. is now a resonance (a pole in the second Riemann sheet) yielding an expo- 
nentially damped amplitude that oscillates with frequency Mq. The damping 
being the width of the resonance [see eq. (|3.53|) 1- 

Eventually, for very large r, the power-like tail coming from the cut contribution 
will dominate over the exponentially damped contribution of the resonance. 

T^O. It is similar as for T = 0, except that the damping rate becomes a non-trivial 
function of t, r, T and M given by eq. (|3.59|) . 

The thermal cut with its non-oscillating contribution to <p(x, t) is also present 
here and dominates for r — > oo. 

All the asymptotic results hold to first order in the field amplitude. However, since the 
field vanishes for r — > oo, our results are true for any initial amplitude provided r is large 
enough. 

The relaxation of an initially Gaussian inhomogeneous fluctuation is described in terms 
of the spreading of the packet and decay in spherical waves, and we found that the time scale 
for relaxation, production and absorption of pions is a function of rapidity such that for larger 
rapidities the relaxational and decay processes are slower. We are currently extending these 
studies to the case of non-equilibrium fluctuations produced during the stage of parametric 
amplification as the a rolls to the ground state |jTo| . 

We believe that the techniques developed in this work and the non-equilibrium aspects 
found here will be of use in a variety of physical contexts. One that comes to mind concerns 
the emission of Goldstone bosons from an axion string |35| , p6| . This emission is important 
in determining the exact upper bound on the axion decay constant, and thus on whether 
axion models for solving the strong CP problem are still viable. This problem will entail 
some modifications of the tools described here, the most notable one having to do with the 
existence of zero modes for the string configuration. 
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Figure Captions 

Fig. 1: One- loop Feynman diagrams contributing to the equation of motion. The dashed 
line corresponds to the insertion of the background field <fi. Thin lines correspond to pion 
propagators, thick lines to a propagators. The tadpoles had been absorbed in mass and f n 
renor malizat ion . 

Fig. 2: Comparison between S(y,p,T) / 'M a and the Breit-Wigner approximation for T = 
O.lMj ; \p\ = 10M a vs y = \Juj 2 — p' 2 /M a . Solid line is the Breit-Wigner approximation, 
dotted line the full spectral density. 
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This figure "figl-l.png" is available in "png" format from: 
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